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Abstract 

This paper describes a technique for analyzing the stochastic structure of the vectorial 
capacity using moment-generating functions. In such formulation, for an infectious disease 
transmitted by a vector, we obtain the generating function for the distribution of the number 
of infectious contacts (e.g., infectious bites from mosquitoes) between vectors and humans 
after a contact to a single infected individual. This approach permits us to generally derive 
the moments of the distribution and, under some conditions, derive the distribution function 
of the vectorial capacity. A stochastic modeling framework is helpful for analyzing the 
dynamics of disease spreading, such as performing sensitivity analysis. 


Introduction 


The field of mathema t ical m odeling of infectious diseases had seminal works by Rossi ( IQld ) 
an d iMacdonald et all (Il957ll that considered t he dy n amics of malaria spreading. The w orks 


by Garrett- JonesI ( 19641 ): Garrett-Jones et ah (1964); Garrett-Jones and Shidrawi ( 19691 ) i 


_ ^ ^ in¬ 
troduced the concept of Garret-Jones’ vectorial capacity, which is the number of future 
inoculations that arise from an infectious case. Such modeling framework has been used not 
only for malaria, but also for other vector-transmitted diseases. This is the case of dengue, 
a disease already well-known for causing epidemics in various parts of the world and also 
emerging infections, such as chikungunya. 


The modeling framework has been described in the literature, for instance bv lAnderson and Mav 


(jl991r ). More recent works also rev i sited t he issue in a tutorial-oriented style or a historical 
perspective. ISmith and McKenziel (120041 ) revisit the formulae involved in the statics and 
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the dynamics of malaria infections. ISmitli et al.l (120121) provide t he history behind the con¬ 
tributions of Ross and MacDonald. iMassad and Coutinhol (j2012l ) provide a treatment that 
explains the dimensionless property of the number. 


Bailey et al. ( 19821 . 19751 ) analyzed epidemics using stochastic models, including using 


moment-generating functions, which provided a seminal body of work in its own right, but 
did not model th e vectorial capacity as p resented here. Also, for a general epidemic pro¬ 
cess, for instance. iBall and Donnellvl (119951 ) presents approximations results, using moment¬ 
generating functions. 


Nevertheless, the treatment presented here using moment-generating functions to de¬ 
scribe the vectorial capacity has not appeared elsewhere, to the best of our knowledge, and 
can be quite important to advance the study of the vectorial capacity. 


Model 


In order to describe the vectorial capacity we need to describe the biting process and the 
survival of vectors. Vectors are species that can transmit pathogens causing diseases. In 
the case of malaria or dengue, vectors are typically mosquitoes. For this reason, we call 
the process between the vector and humans as a biting process, but it could be applied to 
species, other than mosquitoes. 

The total number of vectors is a random variable M and the number of humans is 
assumed fixed at h individuals. 


Moment generating functions 


For any generic random variable X the moment-generating function (MGF) is defined by 
X*{e) = E[e^^]^ Therefore, 


poo 

X*{e)= / e^^P{X = x)dx. 
Jo 


( 1 ) 


For discrete distributions, which is the case of the variable that describes the vectorial 
capacity is convenient to work with the z-transform: 


X*{z) = z^P{X = fc), 2 : = e^, if X is discrete. 


A:=0 


^The Laplace transform is the case when there is a variable replacement from 9 to —s, i.e. E[e . 
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The biting process 


The individual number of bites is described by a random variable B that has distribution 
B[x) and MGF B*{z). When counting the multiple individuals, with the assumption that 
the count numbers are i.i.d, we have a total number Bt = • Hence, the MGF for 

Bt is B^{z) = M{B{z)). 


The survival process 

Let 5 be a random variable that describes the survival time of a female mosquito after biting 
an infected person and the survival time after the incubation period. The distribution 
of S and Sy are given by functions S{t) and respectively. The moment generating 

functions follow from the dehnition S{9) = E[e^^] and Sy{9) = . 


Biting as a Poisson process 


We now introduce an initial assumption that intervals between contacts for any vector are 
independent identically distributed random variables. A second assumption is that the 
distribution of the intervals between contacts is exponential with mean A. 


Let H be a random variable that describes the number of bites by a mosquito in a given 
period of time. Since the intervals between contacts are exponentially distributed, then B 
has a Poisson distribution. Therefore, the number B of bites (individual) over a period t 
has a density distribution given by 


and its MGF given by 


P{B = k) = e 




kl 


B*{z) = 


( 2 ) 


Taking into account all mosquitoes, then the MGF of the total number of bites is the 
product of moment generating function of B expressed in Eq. [2l 

B*^{z) = M{B*{z)) (3) 

where M(z) is the moment generating function of number of mosquitoes. For instance, if 
the total number of mosquitoes is constant at M = m: 

= (4) 
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Now we turn to the number of bites of a mosquito during its virus-carrying state. 
The number of bites in infective state is also described by a Poisson distribution with 
parameter A. If we condition the number of bites to the time S , we have: 

{Bl\S = t} = B ^ Poisson{Xt). 

Following Eq. [2l the moment-generating function is then given by: 

Bl{z\S = t) = (5) 

We find B'^ by integrating = t) over the spectrum of the probability density 

function of S: 

1*00 

Bl{z)= = t)dt = S*{X{z-l)) (6) 

Jo 

Vectorial capacity process 

We start with a single infected human. Let Ns be a random variable that describes the 
number of mosquitoes that bite this infected individual, get blood with virus and also survive 
past the incubation period. Then the total number Bsi of bites given by this group is the 
sum of bites per individual, i.e., Bsi = ^L,j ■ 

Therefore, assuming all B^j are i.i.d. (independent, identically distributed), and using 
Eq. [6l the moment-generating function is 

B*si{z) = N,{S*{X{z-l))) (7) 

Now we have to find the moment generating function of Ns ■ Let the incubation time of 
the pathogen be a random variable V. 

There are a number of bites Bj that effectively infect mosquitoes. Some of these will not 
survive past the incubation time V. Therefore, there is a probability that S >V. From all 
bites that infect mosquitoes, the number Ns is given then by a binomial distribution: 

Ns\V, Bi, S ~ Binom{Bi, P{S > V)) 

where the probability pv that a vector survives past the incubation period is given by 

poo 

PV = P{S >V)= PiS> t)P{V = t)dt. 

Jo 

The moment-generating function for Ns conditioned on Bj , given its binomial distribu¬ 
tion, is 

N*s{z\Bj = k) = {pvz + l-pv)^ (8) 
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The MGF of the number of mosquitoes that reach the state of effective vectors (past the 
incubation period), conditioned oir the total irumber of virus-carrying ones, is given by 

= k) = Ns{S*{\{z - 1))) = {pv{S*{X{z - 1))) + l-pvf. (9) 

Now the number Bj of bites that infect mosquitoes has a binomial distribution, and 
its parameter is the probability q to bite the single human infected. If we assume that 
the probability to bite the single human infected is the same (uniformly distributed) across 
mosquitoes, then q = hi/h, as function of the number hi of infected humans. If hj = 1, we 
have simply q = \/h. 

Bi Binom{BT, q)- 

Therefore, the moment-generating function of Bj can be expressed from the moment¬ 
generating function of Bt derived in Eq. UJ 

B}{z) = B^iqz + l-q)= (iQ) 


Finally, using this MGF of Bj{z) and the result in Eq. EJ 

BMz) = BKpv{S*{X{z-1))) + 1-pv) 

= B^q{pv{S*{X{z-l))) + l-pv) + l-q) 

B*si{z) = M*{B*iq{pv{S*{X{z-l))) + l-pv) + l-q)) (11) 


Eq. [TT] is a more general result to derive the genrating function. In particular, under 
the assumption of a Poisson distribution for the number of bites and a fixed number m of 
mosquitoes: 

Bgj(z) = ( 12 ) 


Lifetime as an exponentially—distributed random variable 

We now assume that lifetime (or at least adult lifetime) is exponentially distributed. Then, 
upon this assumption, S also has an exponential distribution. The probability to survive a 
given period of time, say a day, is p = e~^ and for t days we have p* = 

Therefore, we can express P{S > V) in terms of the moment-generating function of V: 

roo 

pv = P{S>V)= e-<^^P{V = t)dt = V*{e)\e=-g (13) 

JQ 
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Moments of the vectorial process 

If we consider the latent period to be deterministic, given by v days, then V*{z) = 
(discrete case) or V*{9) = (continuous case). 

Also, the lifetime distribution with an exponential distribution has a probability density 
function P{S = t) = ge~^^ and its MGF S*{9) = —gl{9 — g). 

Then we apply S*{9) = —gl{9 — g) and pv = e~^'" into Eq. [121 

B*sj{z) = (14) 


If we take the first derivative of B*gj{z ), for 2 ; = 1, we find the mean E\Bsi] '■ 


dB*si{z), . . m\^e 9^ 

^^U = E[Bal = —^ 

which is the well-known formula for vectorial capacity. 

If we take the second-order derivative, for 2 ; = 1, we find E[Bgj] — E[Bsi] ■ 
(PB%t(z) dB%j m\^e~9'’g mX^e~9'"g 2A 


dz'^ 


h{X{z - 1) -fif)^ 


h {X{z-l)-gY 


d?B%j, ,,2 2mX^e-9'" ,,2 i2A 

2=1 = {E[Bsi])^ H- ^ -= {E[Bgi])^ + E[Bsi] —. 


dz'^ 


hg^ 
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(15) 


Thus, the variance is found: Var{Bsi) = E[Bsi]il + y) = "" (1 + y). 


Inverse transform to obtain the distribution 

Using the power series rewrite Eq. [TTJ 

_ 

k=0 


where f{z) = mX^—^ ^ ■ Rewriting previous equation after some algebraic modi¬ 

fications, 




fc =0 


'-^k\\h{g + X)) 


(16) 


g+\ 
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Table 1: Model: variables and their descriptions 

m Total number of mosquitoes 

h Total number of humans 

A biting rate 

g parameter for the survival function, when exponentially distributed 

V incubation time, when considered deterministic 

B Number of individual bites in a period of time - Poisson distributed 
Bt Total number of bites in the period of time 

Bj Total number of bites that turn mosquitoes into an virus-carrying state 
Bl Number of bites of a mosquito during its period in virus-carrying state 
Bsi Total number of bites by the group of virus-carrying mosquitoes 
Ns Number of survivals of mosquitoes 
S Survival time after biting an infected person 

V Random variable for the incubation time 


Taking the inverse z-transform we find the probability density function for Bsi '■ 


P{Bsi = n) = 



(17) 


where «„ = 1, if n = 0, and = 0, otherwise. This gives us the probability distribution 
for the number of secondary bites that happen when there is a single infected human. 


The power series in the equation converges to aiFi{n + 1,2,a), where a 
and iTi is the confluent hypergeometric function. Hence, 


V h(g+P ) 


P{Bsi =n) = 



X Y / mXe-9^g \ 
g + x) \h{g + X)) 


iFi{n + 1,2, 


/ mXe 

U(5+a)JV 


(18) 


Discussion 


We believe this formulation can be important for studying dynamics of infectious diseases 
transmitted by vectors such as malaria, dengue and chikungunya. 

From our formula for the vectorial capacity process (moment-generating function), we 
are able to obtain all moments of the vectorial capacity process. We demonstrate this in the 
case of simple incubation period and exponential distribution for survival. 

It is desirable to collect data from experiments reported in the literature and possibly 
experiments designed to measure in the field the distributions that describe the biting process 
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and survival process in order to construct generating functions that present good fitting to 
the observed data. We intend to further advance our research in these directions. 
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